%% 用于分析信号的频谱
% 用阴极保护的数据做分析
%% f分析
% 提取数据
load('signal.mat')
% signal=double(S1{1:11300,1});

fs = 1000; % 采样频率，需根据实际系统设置

% 计算 FFT
L = length(signal);
fft_result = fft(signal);
P2 = abs(fft_result / L); % 双侧频谱
P1 = P2(1:L/2+1); % 单侧频谱
P1(2:end-1) = 2 * P1(2:end-1);

% 频率轴
f = fs * (0:(L/2)) / L;

% 绘制频谱
plot(f, P1)
title('FFT of Signal')
xlabel('Frequency (Hz)')
ylabel('P1(f)')


